Self-sustaining dynamical nuclear polarization oscillations in quantum dots 
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Early experiments on spin-blockaded double quantum dots revealed surprising robust, large- 
amplitude current oscillations in the presence of a static (dc) source-drain bias [see e.g. K. Ono, 
S. Tarucha, Phys. Rev. Lett. 92, 256803 (2004)]. Experimental evidence strongly indicates that 
dynamical nuclear polarization plays a central role, but the mechanism has remained a mystery. 
Here we introduce a minimal albeit realistic model of coupled electron and nuclear spin dynamics 
which supports robust self-sustained oscillations. Our mechanism relies on a nuclear-spin analog of 
the tunneling magnetoresistance phenomenon (spin-dependent tunneling rates in the presence of an 
inhomogeneous Overhauser field) and nuclear spin diffusion, which governs dynamics of the spatial 
profile of nuclear polarization. The extremely long oscillation periods (up to hundreds of seconds) 
observed in experiments as well as the differences in phenomenology between vertical and lateral 
quantum dot structures are naturally explained in the proposed framework. 



The coupling of electron and nuclear spin dynamics 
is responsible for a wide variety of intriguing transport 
phenomena in semiconductor devices. Spin exchange 
between electron and nuclear spins provides a mecha- 
nism for electron spin flips which can dramatically alter 
the behavior of systems such as spin-blockaded quantum 
dots, where transport is highly sensitive to spin selec- 
tion rules [THZ]. Furthermore, the nuclear spins produce 
a hyperfine (Overhauser) field that shifts the electronic 
Zeeman energy by an amount corresponding to an effec- 
tive field that may reach as high as a few Tesla when 
the nuclei are fully polarized. This Overhauser field can 
have dramatic consequences for transport in quantum 
dots, where discrete levels may be shifted in-to or out-of 
resonance [H I5HTU] . The combination of these two effects 
- electron-nuclear spin-exchange which polarizes nuclear 
spins, and subsequent back-action on energy-dependent 
spin flip rates - is responsible for a variety of interesting 
nonlinear dynamical effects such as multistability, hys- 
teresis, and intermittencv jSHTS] . 

Among all of the nonlinear phenomena which have 
been observed thus far in transport through double quan- 
tum dots (DQDs), perhaps the most striking is the ap- 
pearance of spontaneous, stable current oscillations un- 
der the application of a dc source-drain biasjll[T^. This 
phenomenon is remarkable for a number of reasons. First, 
the oscillations occur with very long periods ranging from 
seconds to hundreds of seconds. These timescales are 
10'' — 10^ times longer than the (1 pA) /e ~ 100 ns micro- 
scopic timescale associated with single electron tunneling 
through the double dot. Second, the oscillations are ac- 
companied by long transients and long memory times 
when the source-drain bias is switched off and on. Fi- 
nally, after many years of experiments by a variety of 
groups, the oscillations have only ever been seen in ver- 
tical DQDs; the phenomenon has never been observed in 
a gate-defined lateral double quantum dot. 
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FIG. 1: Mechanism of nuclear polarization oscillations in a 
spin-blockaded double quantum dot. Polarization is driven 
on a short time scale by resonant hyperfine transitions inside 
the DQD. Spin injection in the presence of an inhomogeneous 
Overhauser field leads to a polarization overshoot in the dot. 
Nuclear spin diffusion homogenizes the Overhauser field on 
a much longer timescale. Spin-flip transition rates inside the 
DQD adapt but lead to an overshoot in the opposite direction. 

Using nuclear magnetic resonance it was shown that 
the oscillations are in some way driven by nuclear spin 
dynamics[8,. However, despite wide interest in the prob- 
lem, a viable mechanism has thus far remained elusive. 
Here we present a straightforward mechanism which nat- 
urally produces oscillations with similar phenomenology. 
The mechanism relies on nuclear spin diffusion |17|, I18j 
and on spin-dependent tunneling rates [19] which are con- 
trolled by the spatial proflle of the Overhauser field. 

Nuclear spin diffusion, being a slow process, introduces 
the correct timescale into the dynamics. Furthermore, 
it also accounts for a sharp difference in predicted phe- 
nomenology for vertical and lateral DQD structures. The 
length scale for out-of-dot diffusion is set by the combi- 
nation of barrier and quantum well half-widths, and is 
typically a few tens of nanometers. For typical diffusion 
parameters [T71 [T5] this translates into diffusion times on 
the order of 10 seconds, consistent with the observed 
oscillation period values. These timescales are much 
longer than those arising from coherent mechanisms j20| . 
Additionally, in vertical DQDs such as those used in 
Ref. O HI], the edges of the dot are defined by the 
mesa structure itself. In such structures, nuclear po- 
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larization predominantly diffuses in the vertical direc- 
tion, into the adjacent tunnel barriers. For gate-defined 
lateral DQDs, however, nuclear polarization can diffuse 
in all directions [24] ■ and is not expected to flow signifi- 
cantly into the tunnel barriers. Thus the feedback effect 
via polarization-dependent tunneling, which is responsi- 
ble for the oscillations, is expected to occur in vertical 
DQDs but not in lateral DQDs. This is consistent with 
the observation that oscillations are frequently observed 
in vertical structures but never in the lateral structures. 

Schematically, oscillations arise as described in Fig. [l] 
An initial imbalance of hyperfine spin flip rates for up 
and down electron spins leads to a fast build up of nu- 
clear polarization inside the DQD. The resulting inho- 
mogeneity of the Overhauser field between the DQD and 
its surroundings enhances the probability of injecting the 
spins with the dominant hyperfine rate. This causes the 
polarization inside the dot to "overshoot." On a much 
longer time scale, nuclear spin diffusion homogenizes the 
Overhauser field. As the spin-injection probabilities re- 
act accordingly, the balance of hyperfine transition rates 
inside the DQD reverses, and starts to drive the nuclear 
polarization in the dot back toward zero. In a similar 
way, the polarization inside the dot again overshoots and 
then the cycle repeats. 

Here we describe the coupled electron and nuclear 
spin dynamics through the simplest possible model which 
qualitatively captures the behavior of the essential phys- 
ical degrees of freedom. In principle, detailed numerical 
modeling of the dot and of the full spatial profile of nu- 
clear polarization could be attempted. However, such 
modeling would introduce a much higher level of com- 
plexity, and would likely cloud rather than clarify the 
essence of the oscillation mechanism. Instead, we will 
write a set of dynamical equations for two polarization 
variables, one representing the polarization within the 
double quantum dot, and one representing the polariza- 
tion under the tunnel barrier to the source lead. The 
intradot polarization variable is driven by hyperfine spin 
fiip processes with electron spins within the DQD. Polar- 
ization is then transferred to the barrier via spin diffusion 
with a large time constant. The delayed reaction of the 
barrier polarization variable to the intradot spin dynam- 
ics leads to oscillations as outlined above. 

A key to our mechanism is the difference in probabili- 
ties for spin-up and spin-down electrons to tunnel into the 
quantum dot when it is empty [T9]. Naively, one might 
expect the respective tunneling rates to differ due to the 
application of a homogeneous Zeeman field, since the final 
state energies are different. However, as shown in Figj2^, 
up and down spins tunnel under identical Zeeman-shifted 
barriers. Provided that the dot levels are set far below 
the chemical potential of the lead, Ep, and that the lead 
has an approximately constant density of states in the 
energy range of interest, the tunneling-rates for the two 
spin species are equal in the case of a homogeneous field. 



a) Equal Tunneling b) Spin Up Favored 




Vj>K- 



y<0 



c) Spin Down Favored d) Equal Tunneling 







x<0 




¥ 

t" 










Ax 










FIG. 2: Spin-dependent tunneling due to inhomogeneous 
Overhauser field, a) For a homogeneous Zeeman field, up 
and down spins are subjected to identical barriers and tunnel 
into an empty dot with equal probabilities, b) When nuclear 
polarization is nonzero only under the barrier, 1/ 7^ 0, up 
and down spins are subjected to different barriers [B = Q for 
illustration in b-d). c) When nuclear polarization is nonzero 
only inside the quantum dot, a; 7^ 0, up and down spins tunnel 
in at different relative energies, d) If the nuclear polarizations 
in the dot and in the barrier are nonzero but equal, x = y, np 
and down spins tunnel in with equal probabilities. 

What happens in the case of an inhomogeneous 
Zeeman- Overhauser field? For demonstration, consider 
the case shown in Figj2]D, where the nuclear polarization 
is large under the barrier and zero outside. Here the 
Overhauser field locally increases the Zeeman energy un- 
der the barrier, effectively creating a higher barrier for 
down spins, and a lower barrier for up spins. In this sit- 
uation, an empty dot is more likely to be filled by a spin- 
up electron than by a spin-down electron. Similarly, as 
shown in Figj2];, nuclear polarization concentrated only 
inside the dot can also affect the tunneling-in probabil- 
ities by changing the tunneling energies relative to the 
tops of the spin-up and spin-down barriers. 

It is interesting to note the similarity between this ef- 
fect and the phenomenon of tunneling magnetoresistance 
(TMR) [2TI [22] . In both cases transport is dominated by 
tunneling through a barrier and spin polarization is used 
as a knob to control tunneling rate. While in TMR the 
spin polarization is due to magnetization in the regions 
surrounding the barrier, in our case the dominant effect 
is due to under-barrier nuclear spin polarization. The 
discovery of TMR has had important consequences for 
magnetic memory applications. One can envision that 
some of these ideas can be transposed to DQD systems. 

We now consider sequential electron transport through 
a spin-blockaded double quantum dot connected to leads 
with an applied dc source-drain bias, as described for 
example in Refs. [SI [11] [T^]. In the two-electron spin- 
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FIG. 3: Two-electron energy levels involved in spin- 
blockaded transport (adapted from Ref . [25] ) . a) As a func- 
tion of interdot potential bias, which controls the asymmetry 
of the double well potential, the (1,1) and (0,2) singlet states 
exhibit an anticrossing. b) Energy levels at large detuning, 
indicated by the dashed vertical line in a). The singlet lev- 
els are broadened due to the coupling of the (0,2) state to 
the drain lead. Hyperfine-assisted transitions from |r±) to 
IS) provide a source for the nuclear polarization x within the 
double dot. 

blockade regime, "(1,1)" orbital configurations with one 
electron in each dot, and a "(0,2)" configuration with 
both electrons in the second dot have nearly the same 
electrostatic energies. In the (1,1) configuration, where 
overlap between electrons is negligible, all four spin states 
(one singlet and three triplet states) are nearly degener- 
ate in energy. For the (0,2) configuration, however, only 
the spin singlet configuration is allowed due to the Pauli 
exclusion principle (the single dot orbital level spacing is 
assumed to be much larger than the applied bias) . Inter- 
dot tunneling hybridizes the (1,1) and (0,2) singlet states, 
producing the states labeled \S) and \S') in Fig. |3] 

Tunneling out of the double dot occurs from the (0,2) 
singlet state, which is coupled to the drain lead. Through 
hybridization, both singlet states \S) and |S") acquire 
finite lifetimes, reflected in their broadened lineshapes as 
shown in Fig. [3]3. When only spin-conserving tunneling 
processes are taken into account, the triplet states remain 
decoupled from the drain. Therefore the rate-limiting 
step which controls current through this system is the 
decay of the long-lived triplet states through resonant 
hyperfine-assisted transitions to the singlet states jiS*) and 
1 5'), or higher order processes which may also break the 
conservation of spin within the double dot 6 . Hyperfine 
assisted transitions from \T±) to \S) and \S') transfer 
angular momentum from electron to nuclear spins, and 
thus drive the nuclear polarization dynamics. 

Here we focus on the regime of large detuning where 
the level \S') is far separated in energy from the triplet 
states and can be ignored in the calculation of hyperfine- 
assisted triplet-singlet transitions. We seek a coupled 
set of dynamical equations in two polarization variables. 
The first variable, x = {N+ — N-)/{N++N-), represents 
the fractional polarization inside the double dot, where 
N+ {N-) is the number of nuclear spins oriented along 
(against) the external field. For a typical device, N = 



N+ + N_ 



10 . The second variable, which we denote 



by y, represents the fractional polarization within the 
tunnel barrier connecting the source lead to the first dot. 



The intradot polarization x controls feedback through 
the Overhauser shift of the electronic triplet levels, which 
can bring these levels into or out-of resonance with the 
singlet. The energies e± of the triplet states |T±), relative 
to 1 5), are given by 



e± = e± g*fiBB ± Ax, 



(1) 



where e is the singlet-triplet detuning, g* is the electronic 
effective g-factor {g* « —0.44 in GaAs), /ib is the Bohr 
magneton, B is the strength of the applied magnetic field, 
and A ~ 100 fj,eV is the hyperfine coupling constant. 

Each time an electron decays from |T+) or |r_) to \S) 
via hyperfine exchange, one nuclear spin is flipped from 
down to up, or up to down, respectively. The probability 
for an electron that enters the dot to cause a positive 
(negative) increment to the nuclear polarization during 
its escape is determined by the probability /+ (/-) that 
the electron entered into the state |T+) {\T_)), and by 
the probability that the electron escapes via the hyperfine 
exchange process rather than by alternative nuclear-spin- 
independent mechanisms [111 I23| . The hyperfine spin flip 
probabilities are determined by the ratios Vl^]^^/(Vl^|^^ + 
VF'"), where W^^ is the hyperfine decay rate of jT-t) and 
VF'" describes the collective effects of spin-orbit coupling, 
spin exchange with the leads, and cotunneling processes. 

In our model, we assume that all nuclear spin flips due 
to hyperfine exchange with the electrons occur within the 
double quantum dot. Therefore the dot polarization x 
receives kicks (with magnitude 1/N) on the timescale 
of single electron hopping through the dot, 100 ns to 
1 /is, while the barrier polarization y has no dynamics 
on this small timescale. On a much longer timescale, 
nuclear polarization may diffuse from the dot region into 
the barrier region, providing a source for y. 

Mathematically, it is simplest to analyze the regime 
where ^ Wl'^. Here the total current, i.e. the ef- 
fective frequency of electrons passing through the dou- 
ble dot, is determined by W^'". Additionally, the hyper- 
fine decay probabilities reduce to W^^/iy". The de- 
pendence on M^™ cancels from the nuclear polarization 
rate, which depends on products of attempt frequencies 
and spin flip probabilities, leaving behind contributions 
proportional to the hyperfine rates W^^ weighted by the 
loading probabilities f± : 

y = -2TDy 



- /_I^"F 
Tdx, 



(2) 
(3) 



where ~ 0.1 s^^ is the inverse of the time constant 
for diffusion from the dot to the barrier. The hyperfine 
spin-fiip rates W^^ are given by Fermi's golden rule|ll): 



A^{1_ 



-x)j 

+ j2 ' 



(4) 



where 7 is the decay rate of \S) due to its coupling to 
the drain. We account for the dependence of the loading 
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probabilities f± on the Overhauser field inhomogeneity 
in a lowest-order expansion in x and y: 



f± = \[i±vix~y)i 



(5) 



where ry controls the sensitivity of the loading probabili- 
ties to a polarization gradient. The factors of 2 in front 
of To in Eqs. ([2| and ([s]) account for the fact that polar- 
ization diffuses in both directions (up and down). 

Under what conditions might we expect to find oscilla- 
tions in the flow defined by Eqs. ^ and (|3|? Typically, 
oscillations are found when the linearized system has the 
form of an "unstable spiral:" 



(6) 



with (a + /3) > and (a — /3)^ — 4^ < 0. These conditions 
ensure that the eigenvalues are complex, with positive 
real part. Comparing with Eq.([3]), we see that y ~ x, 
with a positive coefficient of x due to the fact that po- 
larization preserves its sign as it flows into the barrier. 
Therefore, we need the coefficient of y in Eq.([2]) to be 
negative. Substituting expression ([s]) for f± into Eq.Q, 
we thus obtain a threshold rj > ATd/Wq, where Wq is 
the hyperfine spin flip rate at the unstable fixed point. 

Going one step further, we can expand Eqs. ^ and 
([3]) in the deviations x and y from the (unstable) fixed 
point of the nonlinear system. Notably, because y only 
appears to linear order in the original expressions, only 
y-independent or y-linear terms will show up in the ex- 
pansion. In general, all other terms will appear: 



X w ciqx + coiy + 



y = Tdx -~ 2TDy. (7) 



where the dots represent higher order terms C2ox^ + 
ciixy + C3oX^ + C2ix^y+- ■ ■ . Comparing to Eq.(|6]), we see 
that we need cio > 2Td > to ensure a positive real part 
of the eigenvalues, and cqi < 0, (cio -t- ^Td)^ < ^coi\Td 
to ensure a negative discriminant. These considerations 
lead to the oscillatory regime shown in Fig. |4] 

Using the vast separation of timescales between the 
hyperfine spin-flip driven polarization dynamics and the 
slow diffusion processes, we explore another avenue of 
analysis. We assume that the barrier polarization y is 
constant on the timescale of changes in the dot polariza- 
tion and examine the flxed points of the resulting quasi- 
one-dimensional dynamical system ([2| . Figure |4] shows 
the corresponding stable (green) and unstable (red) flxed 
points, superimposed on the velocity field map of the full 
system (arrows indicate direction of the polarization ve- 
locity, and the color scale indicates its magnitude). The 
oscillations can be seen as resulting from the combination 
of the circulatory flow around the origin, ensured by the 
spiral condition above, combined with the existence of 
the unstable branch of quasi-fixed-points near the origin. 
Fast horizontal motion towards the quasi-stable points, 
followed by slow drift along the green curves provides a 
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FIG. 4; Polarization velocity field in the oscillatory regime. 
Parameter values: e/A ^ -3, 7M = 0.05, 5 = 0,?) = 0.4, 
Td/A = 10~^^. Arrows indicate the direction of the velocity 
field (i, y), while the color scale indicates the magnitude of the 
velocity (arbitrary units). The green and red curves indicate 
branches of stable and unstable fixed points of the quasi-one- 
dimensional dynamics ([2| with y held constant. The blue 
curve shows the approximate trajectory of the limit cycle. 

pictorial representation of the oscillation mechanism out- 
lined in Fig.[l] The oscillation period is dominated by the 
length of the excursions along the green branches. As a 
result, the period grows as the oscillation amplitude in- 
creases, consistent with experiment [5]. As parameters 
vary, a transition out of the oscillation regime occurs 
when the unstable branch shrinks and disappears. 

We have identified a straightforward physical mecha- 
nism which can produce stable oscillations of dynami- 
cal nuclear polarization in spin-blockaded vertical double 
quantum dots. The mechanism relies on nuclear spin 
diffusion into a tunnel barrier and is active only in ver- 
tical DQDs. The dependence of spin-injection probabili- 
ties and spin diffusion times on barrier width provides a 
clear experimental signature of this mechanism. Persis- 
tent oscillations can serve as a new probe of nuclear spin 
diffusion and spin dynamics in vertical structures. 
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